Gravitational collapse in braneworld models with curvature corrections 
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We study the collapse of a homogeneous braneworld dust cloud in the context of the various 
curvature correction scenarios, namely, the induced-gravity, the Gauss-Bonnet, and the combined 
induced- gravity and Gauss-Bonnet. In accordance to the Randall-Sundrum model, and contrary to 
four-dimensional general relativity, we show in all cases that the exterior spacetime on the brane is 
non- static. 
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In this Letter, we discuss the Oppenheimer-Snyder-like 
collapse on a brane in the context of curvature correction 
terms. In the Randall-Sundrum scenario, this problem 
has been analyzed in 1], and found that, contrary to the 
general relativity case, the vacuum exterior of a spheri- 
cal cloud is non-static. This is a result of modification of 
the effective Einstein equations on the brane with local 
and non-local terms representing high energy corrections 
to general relativity. The non-static nature of the exte- 
rior metric mainly arises because of the presence of bulk 
graviton stresses, which transmit effects non-locally from 
the interior to the exterior on the brane, and of the non- 
vanishing of the effective pressure at the boundary sur- 
face which connects the interior with exterior metric 
via the four dimensional matching conditions. 

We derive here the same result within the induced- 
gravity, the Gauss-Bonnet, and the combined Gauss- 
Bonnet and induced gravity braneworld models. In all 
these models, the effective Einstein equations on the 
brane are modified by local and non-local terms, which 
are much more complicated than the corresponding terms 
in the Randall-Sundrum case, but nevertheless, the non- 
staticity arises because of a mismatch of the interior with 
the exterior metric on the boundary collapsing surface. 

We consider for convenience and without loss of gen- 
erality the extra-dimensional coordinate y such that the 
brane is fixed at y = 0. The induced metric h^^, on 
this hypersurface is defined by Hab = 9ab — nAns, 
with the unit vector normal to the brane (^, v = 
0, 1, 2, 3; A, B = 0, 1, 2, 3, 5). The total action of the sys- 
tem is taken to be 



S=^ J d'^xy^ I 7^ - 2A, 



5 + a 



} 



where TZ, R are the Ricci scalars of the metrics gAB and 
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Hab respectively. The Gauss-Bonnet coupling a has di- 
mensions (length)'^ and is defined as 



1 



(2) 



with gs the string energy scale, while the induced-gravity 
crossover lenght scale Tc is 



Mi 
Mi 



(3) 



Here, the fundamental (M5) and the four-dimensional 
{M4) Planck masses are given by 



kI = SttGs = M^^ , kI = 87rG4 = M4 
The brane tension is given by 



A4 



(4) 



(5) 



and is non-negative. (Note that A4 is not the same as 
the cosmological constant on the brane.) 

The collapse region has in comoving coordinates a 
Robertson- Walker metric 

= -dr^ + a{rf{l + kxV {dx' + x'dill), (6) 

where the scale factor a(T) is given by the modified Fried- 
mann equation of the corresponding model, while the en- 
ergy density is given by the usual dust law p — po{aQ / a)^ , 
with oq standing for the epoch when the cloud started to 
collapse. This Friedmann equation can also be written in 
terms of the proper radius from the center of the cloud 
r(T) = a(T)x/(l-|-fcx^/4) of the collapsing boundary sur- 
face at X = Xo- 

Concerning the exterior of the collapse region, the most 
general static spherically symmetric metric is written in 
standard coordinates as 



ds^ = -F{rfA{r)dt^ + Airy^dr^ + r^dn 



2 • 



(7) 



In order a metric of the form {T)) to be the exterior of the 
interior metric ^ , the metric and the extrinsic curvature 
have to be continuous across the collapsing boundary sur- 
face. Following the method appearing in we first write 



2 



the standard radial geodesic motion of the freely falling 
boundary surface for the exterior metric 



= -A{r) 



E 



F{rY ' 



(8) 



where the dot denotes derivative with respect to proper 
time T, and i? is a constant. Secondly, transforming to 
null coordinates (w,r), where dv = dt + dr/[F(r)A{r)], 
the exterior metric Q takes the form 

ds^ ^ -F{rfA{r)dv^ + 2F{r)dvdr + r^dfll , (9) 

while the interior metric becomes 



ds = — T dv - 



2ar , 



Comparing Eqs. 0, we obtain 



(10) 



(11) 



where E = -fcxo/(l + fcxo/4)^. Then, from Eqs. iP, 
we obtain that F{r) is a constant, and by choosing 
E — 1 + E, we take F{r) = 1. Finally, the candidate 
exterior metric (UJ becomes 

ds^ = ~A{r)dt^+A{r)-^dr^+r^{d0^ + sin^9d<l)^), (12) 

where A(r) is given by Eq. (|lll) . with provided by the 
Friedmann equation of the interior region. 



A. Induced gravity 



The scale factor air) is given by the modified Fried- 
HlllEllllllla of induced gravity 



mann e quat ion 
(a^O) 113 



0" 



± 



AkKp + X) - 2A5 + 



12 12c 



, (13) 



where C is integration constant related to the mass of the 
bulk black hole. We can write Eq. H13|) in terms of the 
proper radius r^r) as 



K4A 



kIiti 1 



E 



±- 



(2kIX - A5 



6X4 AkItti 

Ir^H = — r-4q 



(14) 



where m = PoaoXo/(l + fcXo/4)^ q = Cx^/(1 + kxl/4)\ 
Thus, from Eq. Ijlll) . we obtain 

1 / 



A(r) 



1 

T — 



1 



3 r 



-(2«:2A-A5 + ^)r4 + ^r-4, 



(15) 



For the Randall-Sundrum model, the generic four- 
dimensional effective equations were derived in [lll | . For 
the induced gravity model, such effective braneworld 
equations were derived in (12| in the form 



= kI TI; ~ (k^ A + ^) <5M + ^ (^T^M + ^ St;) , (16) 

where the quantities are given by the algebraic equa- 
tion 



L 



dvdr + r'^dfll (L = L[l)with 



= («2^--A5)<5^;-^2j._^M_^M. (18) 

is the braneworld matter content, if any, while the 
electric part S^^ = C^^gU^n^ of the 5-dimensional Weyl 
tensor C^^^^ carries the influence of non-local gravita- 
tional degrees of freedom in the bulk onto the brane, 
making the brane equations (|16|) not to be, in general, 
closed 13j. Additionally, the energy- momentum tensor 
was shown to satisfy the usual conservation equations 



(19) 



and thus, the Bianchi identities on the brane give from 
Eq. differential equations among : 



= 0, 



(20) 



(semicolon means covariant differentiation with respect 
to h^,^). 

For spherically symmetric braneworld metrics of the 
form ifT^ . the system of Eqs. was fully in- 

tegrated in vacuum in |l4j . [TSi] . and the results are 
as follows. For £^ = on the brane, the solution is 
Schwarzschild- {A)dSi 



A{r) 



1-1 



(21) 



where 7 is integration constant and a = k|A/3 + 2/r'^ — 
2y2^^|A^^A7T67^/%/6rc. For 5^ ^ 0, there are two 
classes of solutions, given in parametric form. The first 
one (with constant is 

A = 1 - 2 - 
r 

5 r 128 ^ /15 23 \ 

+ sqiO- iFi — , — \sq{C)z]z 

r L 105 8 ' 8 ' J 



zs, (22) 
(23) 



3 



where 5 > 0, 7 are integration constants, and a = 
KlX/3 + 2/rl , C = 8(2K|A-A5 + 6/r2)/9r2. The second 
solution (with non-constant Sjf) is 



(24) 



± ^/|.-V3|-^(. + ^/3r 



3(3+^3) V dv 



\v-3\y^ ' 



( S ^k |. - ^/3| 

IW-3IV4 |^, + ^|(V3-l)/8 ' 



where 5 > 0, 7 are integration constants, and a = 
Kj\/3 + 2/rl C = 9(2k|A - A5 + Q/rl)/rl. [The ± sign 
of Eq. (122) is independent of that in Eq. {THJl] . Note that 
the sohitions (|21|I - H25|I are the generic braneworld solu- 
tions and have been obtained without any assumption for 
the bulk space. It is now obvious that the only possible 
static exterior solution ifT^ . surrounding the collapsing 
region, cannot take one of the permissible forms (|21|l . 
(|22|l . H24(l of induced- gravity theory, which means that 
the no-go theorem for induced gravity has been proved. 

We are going to give now another way of showing 
the no-go theorem of induced gravity, without using the 
above static sohitions of the model. This is necessary in 
the case that such exact solutions are not known, as e.g. 
in the Gauss-Bonnet or in the combined Gauss-Bonnet 
and induced gravity braneworld. For the induced-gravity 
scenario, alternatively to Eqs. H16() . generic covariant ef- 
fective braneworld equations have been derived in in 
the form (for the vacuum case) 

(l^^r,4)Gi ^ ~\{h, + ^ ni) + rliri - (26) 



where 
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It is obvious that by taking the trace of Eqs. (|26|) . the 
electric components of the Weyl tensor, which from the 
braneworld viewpoint are not intrinsic quantities, disap- 
pear, and the resulting equation is 



B. Gauss-Bonnet term 



The cosmology of the braneworld Gauss-Bonnet model 
(r^ 0) is given by the Friedmann equation p^ll^ 



/a\^ k 1 / „ (347" 
(a) ='^+8^ -2 + ^ + ^ 



64/^ 



(28) 



where the dimensionless quantities /, J are given by 



8 



, 4 , 8aC 
1 + -aAs + — 
3 a* 



1/2 



(29) 



J 



V2 V 2 

Similarly to the induced gravity case, we obtain 

64/(r)2 



2/3 



.(30) 



where 
/(r) = 

J(r) = 



A{r) = 1 + 



4 . 
1 + g^As 



8a 



J(r) 



J(r) 



8aq 



1/2 



(31) 
(32) 

2 
■3 

(33) 



The generic effective braneworld equations for the 
Gauss-Bonnet model have been derived in 18] in the form 

a(10A5-«I)/ _ 1 J ^ A5 
2(3 -I- aM) V 4 ^ / 4 ' ^ ' 

where 

-i V(7M2-24A/„/3Af"^ + 3A/„^^5Af"'5'^^) , (35) 



- + ^ (1 + ^^-5) - ^ (a5 + ^.^) = . 

(27) 

We can now check by substituting the candidate solution 
(|15|l into Eq. H27I) that this is not satisfied, which means 
that the interior solution cannot match to any static ex- 
terior. 

Note that, while the Friedmann equation H13|l was ini- 
tially derived in Q under the assumption of a particular 
bulk ansatz inducing, of course, the metric ((HJ on the 
brane, however, the derivation of this Friedmann equa- 
tion was also given in in a bulk-independent way, thus, 
showing the full generality of the cosmology p3|l . 



+ ^h^,Mp,EP'' + 3ME^,, (36) 

(3) ^ -AN,,N,+ANP{Np^, + Np,^) + 2Np,^NP% 
+mf,p,Nf^+3h^,{Na,N"- liV^^^TV"^''), (37) 

I = - SMapM"''^ + MapjsM"''^'^^ - 8NpNP 

+mp,,NP'''' - UMp^EP^ . (38) 
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These equations contain the quantities Mafs-ySj ^fiup, 
which are expressed in terms of the induced metric h^^, 
and the extrinsic curvature Kfj^^, as fohows 



Mafias — Rafi-fS — Ka^Kps + KasKfj^ 



(39) 



M, 



a/3 



07/3(5 



M 



h'^^M^p (40) 



up: 



K^p,, , Np^hP'^Np^,, (41) 



where Kp^, satisfies the matching conditions J^J (for the 
vacuum case) 



2n/ 



S':GP)K^p] 



nix 



5^ 



(42) 



with 



ppp ^ j^pp 



25P R" 



(J = J^). It is remarkable that despite the complexity 
of the above equations, by taking the trace of Eq. H34|) 
the quantities Ep^, disappear (as well as the quantities 
Npi,p containing covariant derivatives of Kp^, with re- 
spect to hp^), and the resulting equation is purely four- 
dimensional 

M + a{M'^ - iMapM"^ + Maff-ysM"'^^^) = 2A5 . (43) 

After solving the cubic system H42(l for Kj^, one has 
Eq. constructed solely out of the induced metric hp^. 

In order to do so, it is convenient to write the metric 
(|12|l in a form where its angular part appears as a two- 
dimensional conformally Euclidean space 



-A{r)df 



+ f{xi,X2){dxl+dxl), (44) 



where f{xi,X2) = [1 + {xf + a;2)/4]^^, and xi — 
2 tan(0/2) sin 0, X2 = 2tan(6'/2) cos0. Due to the sym- 
metry of the metric (I44|l . the components of Kj^ in the 
coordinates {t,r, xi, X2) take the form 



Kl^ ^diag{Ki,K2,K3,K3 



(45) 



Subtracting, now, the tt-equation from the rr-equation 
of the system H42f) . we obtain the separable equation 



[4:ar^K^ + Aa{A - 1) - r^] {K2 - isTi) , 
with solutions 



± 



1 



64/(r) 



J{r) 



J{r) 



1/2 



K2 = Ki. 



(46) 



(47) 



(48) 



For the solution l(17|l , plugging back into the system iB^ , 
we obtain a system of two equations for Ki , K2 , from 
where the following equation arises 

2aA' - r T V^KlXr^ /2^r'^ + 4a(l - ^) = , (49) 

(prime means differentiation with respect to r). This 
equation is easily seen to be inconsistent with A(r) given 
by Eq. H31|) . For the solution (|48|l . plugging back into 
the system (|42f) . the situation is more complicated, and 
a polynomial equation of fifth degree on Ki arises: 



K3 



4(1 - 2ar-M')i^i -f K5A 



(50) 



SaKf + 2{2aA" - 1) ' 
ar^Kl + a^Kl + a^Kf + a2K^ + aiKi + oq = 0,(51) 

where 

05 = 2q:P + 4q!(1 - A)] , a4 = aKgAr^ 
03 = {2aA" - l)p + 4q!(1 - A)] 
a2 = KlXr[aA' + r{aA" - 1)] 
2ai = 1 - {2a A" - l)[ArA' + A{2aA" ~ 1) - iaA'^] 

+A"[r^ - 4a + a{r^ + 4:a)A"] - {3 + an^X^y^ /Aa 
ao = KlXr{2aA" - l){2arA" +r- 4aA')/16a. (52) 

Supposed that equation (|43|) is valid, for the metric 
(I44II . and after substituting the values of K2,K3 from 
Eqs. (|^ . l|Sn|) . it becomes a polynomial equation of sixth 
degree on Ki 

beKf + biKf + bsKf + b2Kf + bo = 0, (53) 

where 

66 = 32a'^[r'^ +Aa{l ~ A)] 
64 = 16a[(2aA5 - 3^ - 6a - 12aA'{aA' - r) 
+3a{r^ + 4:a)A" - 6aA{2aA" - 1)] 
63 = 32aKlXr{2aA' - r) 

y = 12a-(3 + 3aK^A2 + 8aA5)r2-12aA(2aA" - 1)^ 

+Aa^A"[A{A5r^ - 3) + 3{r^ + 4a) A"] 
60 = {2a A" - 1 ) {2a(r2 + Aa)A"^ - A" [8a + (1 - 4aA5)r2] 
-2(2ayl"-l)[A(2aA"-l) + 2yl'(aA'--r)] 
-|-2-(2A5 + .jtAV2)r2}. (54) 

Equations |(SIJ, have to be satisfied simultaneously. 
After some algebraic manipulations, this system of equa- 
tions is written equivalently as the following system 



F2Kf+FiKi + 1 = 0, 



C2 — (Ci — Fi)Fi — F2 
(Ci — Fi)F2 — C3 



(55) 
(56) 
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where 

(Ci , C2 , C3 ) = (B2 (Si -pi) +|J3 - ^3 , S3 (Si -pi) +P4 - B4 , 

S4(Si-pi)+p5)/(Sl(Sl-pi)+P2-S2), 

(Fi , i^2 ) = ((Cl C2 - C3) (Si -Pi ) - C2 (S2 -P2) + S4 -P4, 

CiC3(Si-pi)-C3(S2-p2)-p5)/ 

((C2-C2)(Si-pi)-Ci(S2-p2)+S3-p3), 
(Si , S2 , S3 , S4) = (piP2-P3 + '73 , PlP3-P4 + qA , 

piPA-Pb, piPd+qe) /{pl-P2+q2), 

and the various pi, qi are related to a^, bi as 

p^ = ai/ao , q^ = bi/bo. (57) 

It is now straightforward to check that the value Ki of 
Eq. (|^ does not satisfy Eq. which means that the 
no-go theorem of the Gauss-Bonnet braneworld has been 
proved. 

C. Gauss-Bonnet and Induced gravity 

The cosmology of the combined Gauss-Bonnet and 
induced grav ity braneworld is given by the Friedmann 
equation |20| 



(- 



4-3/3 



— VS2^6Qcos e±- 
Spa V 3 



(58) 



where the dimensionless quantities P, Q, Q are given 

by 

(59) 



P 

Q 

e{p,Q) 



9r2 ' 
1 + 3/3/ , 



= 13 



1 



4+^+^(P + A) 



- arccos 
3 



2P3 + 27Q2 _ I8PQ 



(60) 
(61) 

(62) 



2(P2 _ 6Q)3/2 

The ± sign in Eq. H58|l is the same as that in Eq. (fT^ . 
The region in (P, Q)-space for which Eq. is defined, 
is 

1 < < ^ , (63) 
2[ 9P - 8 - (4 - 3P)^/2 ] < 27Q 

< 3P[3- V3(3-2P)]. (64) 

From the above equations, we can write the candidate 
black hole metric as 



A{r) = 1 - 



4-3/3 
12/3a ' 



2r2 
3/3a 



VP(r)2 - 6g(r) cos (e{r) ± ^) , (65) 



where 
P(r) = 1 + 3/3/(r) , 



g(r) -/3 



9(r) = - arccos 



■2P(r)-^ + 27Q(r)2 - 18P(r)Q(r) 



2(P(r)2 -6Q(r))3/2 



(66) 
(67) 

.(68) 



Similarly to the induced gravity case, in the pure Gauss- 
Bonnet model, the Friedmann equation H28|) was ob- 
tained in |18| in a bulk-independent way, showing, 
thus, that this is the most general Robertson- Walker 
braneworld cosmology of the model. For the combined 
Gauss-Bonnet and induced gravity braneworld scenario, 
its Friedmann equation H58I) has not been derived yet in 
a bulk-independent way, but we have no reason to expect 
any discrepancy in this case too. 

In the present case, Eqs. (jSHl-lEIl) - as weh as Eq. (gSJ 
- of the pure Gauss-Bonnet case remain unchanged, since 
they are described by bulk information, while the match- 
ing condition 142|) is now modified by setting its right 
hand side equal to 



6 



3r, 



R 
6" 



Si: 



(69) 



Following the same steps as before, we find the same 
separable equation (|46|l . with solutions now 



K. = ±- 



1 - 2 v/S(r)2 - 6Q(r) cos (o (r) ± ^) ] ',(70) 



K2 = Ki. 
Solution H70|l leads to the equation 



(71) 



2a A' -rTVa[i4Xr^ - ( 1 - A-rA')]/2^yr^ +Aa{l~ A) 

= 0, (72) 

which is inconsistent with Eq. (|65|l . Solution H71|l leads 
to two equations of fifth and sixth degree on Ki of the 
form (|51|l and (|53|l respectively, with the corresponding 
coefficients defined now as follows (denoted with primes) 

flg — 05 , 04 = 04 — arc{l — A — rA') , Og = 03 
4 = 02 + ra[l + 2aA'^ + A{2aA" - 1) + rA'{3aA" - 2) 

-A"{2a + ry2)]/2 
a[ =ai- rc{2A' + rA")[4K^Ar + rc{2A' + rA")]/32 
a'n = ao 

b'e = b6 , 64 64 

63 = 63 - 16arc{r - 2aA'){2A' + rA") 

b'^ = b2- 3arc(2A' + rA")[4K^Ar + rr{2A' + rA")]/2 

b'o = bo- (73) 

Equations remain the same for the above 

primed quantities a'^,b'^, and thus, their incompatibility 



6 



can be easily checked, proving the no-go theorem of the 
combined Gauss-Bonnet and induced gravity braneworld. 

Our analysis of all the above considered models is 
based on 4-dimensional solutions or 4-dimensional ef- 
fective braneworld equations, and we have not studied 
the bulk extension of the considered braneworld regions 
(static exterior and Robertson- Walker interior). We 
know however, that for a given continuous boundary met- 
ric and extrinsic curvature, the propagation of the field 
equations in five dimensions is a well defined initial value 
problem, solvable in principle. 

In conclusion, we have studied the Oppenheimer- 
Snyder-like collapse on braneworld models with curva- 
ture corrections. In all cases considered, using the four- 
dimensional effective equations, and without making as- 
sumptions about the bulk, we have found that the ex- 
terior vacuum spacetime on the brane is non-static. We 



have not found the exterior metric, thus, we are not in 
position to know if the gravitational collapse on the brane 
leaves at late times a signature in the exterior, or if, on 
the contrary, the non-static exterior is transient, tending 
to a static geometry. 



Acknowlegements 

We wish to thank R. Maartens for useful discussions. 
G.K. acknowledges partial support from FONDECYT 
grant 3020031 and from Empresas CMPC. Centro de Es- 
tudios Cientificos is a Millennium Science Institute and 
is funded in part by grants from Fundacion Andes and 
the Tinker Foundation. E.P. is partially supported by 
the NTUA research program "Thalis" . 



M. Bruni, C. Germani and R. Maartens, Phys. Rev. Lett. [16] 
87 (2001) 231302 gr-qc/0108013 . 

M. Govender and N. Dadhich, Phys. Lett. B538 (2002) 
233 hep-th/0109086 . 

C. Deffayet, Phys. Lett. B502 (2001) 199 [17] 
hep-th/0010186 . 

H. Collins and B. Holdom, Phys. Rev. D62 (2000) 105009 
hep-ph/0003 1 7 3). 
Y. Shtanov, hep-th/0005193 

N.J. Kim, H.W. Lee and Y.S. Myung, Phys. Lett. B504 

(2001) 323 hep-th/0101091 . 

E. Kiritsis, N. Tetradis and T.N. Tomaras, JHEP 03 

(2002) 019 hep-th/0202037 . 

V. Sahni and Y. Shtanov, Int. J. Mod. Phys. 11 

(2002) 1 gr-qc/0205111 ; U. Alam and V. Sahni, 
astro- ph/0209443 

K. Maeda, S. Mizuno and T. Torii, Phys. Rev. D68 

(2003) 024033 gr-qc/0303039 . 

G. Dvali, G. Gabadadze and M. Porati, Phys. [18] 
Lett. B485 (2000) 208 hep-th/0005016 ; G. Dvali 
and G. Gabadadze, Phys. Rev. D63 (2001) 065007 [19] 
hep-th/0008054 . 

T. Shiromizu, K. Maeda and M. Sasaki, Phys. Rev. D62 
(2000) 024012 gr-qc/9910076 . 

G. Kofinas, JHEP 0108 (2001) 034 hep-th/0108013 . 
R. Maartens, Phys. Rev. D62 (2000) 084023 
hep-th/0004166 ; R. Maartens, gr-qc/0101059 
G. Kofinas, E. Papantonopoulos and I. Pappa, Phys. Rev. 
D66 (2002) 104014 hep-th/0112019 . [20] 
G. Kofinas, E. Papantonopoulos and V. Zamarias, Phys. 
Rev. D66 (2002) 104028 hep-th/0208207. . 



C. Charmousis and J. Dufaux, Class. Quantum Grav. 
19 (2002) 4671 hep-th/0202107 ; C. Germani and 
C . Sopuer ta, Phys. Rev. Lett. 88 (2002) 231101 
hep-th/0202060 . 

N. Deruelle and T. Dolezel, Phys. Rev. D62 (2000) 
103502 gr-qc/000402j]; I. Lo w and A. Zee, Nucl. 
Phys. B585 (2000) 395 |hep-th/0004124 ; LP. Neupane, 
JHEP 0009 (2000) 040lhep-th/0008190 ; B. Abdesse- 
1am and N. Mohammedi, Phys. Rev. D65 (2002) 084018 

hep-th/0110143 ; Y.M. Cho and LP. Neupane, Phys. 
Rev. D66 (2002) 024044 hep-th/0202140 ; J.E. Lid- 
sey, S. Nojiri and S. Odintsov, JHEP 06 (2002) 026 
[hep-th/0202198 ; J.E. Lidsey and N.J. Nunes, Phys. 
Rev. D67 (2003) 103510 astro-ph/0303168 ; J.P. Gre- 
gory an d A. Padilla, Class. Quant. Grav. 20 (2003) 4221 
"hep-th/0304250 ; C. Barcelo, C. Germani and CP. Sop- 
uerta, Phys. Rev. D68 (2003) 104007 gr-qc/0306072 ; C. 
Charmousis and J. Dufaux, hep-th/0311267 
K. Maeda and T. Torii, Phys. Rev. D69 (2004) 024002 

hep-th/0309152 . 

S. Davis, Phys. Rev. D67 (2003) 024030 
hep-th/0208205 ; E. Gravannis and S. Willison, 

Phys. Lett. 562B (2003) 118 hep-th/0209076 ; N.E. 

Mavromatos and J. Rizos, Phys. Rev. D62 (2000) 124004 
hep-th/0008074 ; Y.M. Cho, LP. Neupane and P.S. 

Wesson, Nucl. Phys. B621 (2002) 388 hep-th/0104227'; 

S. Nojiri and S.D. Odintsov, JHEP 07 (2000) 049 

"hep-th/0006232 . 

G. Kofinas, R. Maartens and E. Papantonopoulos, JHEP 



10 (2003) 066 |hep-th/0307138| . 



